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Abstract 

We are interested in strong approximations of one-dimensional SDEs which 
O f have non-Lipschitz coefficients and which take values in a domain. Under a 

D ■ set of general assumptions we derive an implicit scheme that preserves the 

^0 ' domain of the SDEs and is strongly convergent with rate one. Moreover, 

we show that this general result can be applied to many SDEs we encounter 
in mathematical finance and bio-mathematics. We will demonstrate flexibil- 
ity of our approach by analysing classical examples of SDEs with sublinear 
coefficients (CIR, CEV models and Wright-Fisher diffusion) and also with 
superlinear coefficients (3/2-volatility Ait-Sahalia model). Our goal is to jus- 
tify an efficient Multi-Level Monte Carlo (MLMC) method for a rich family 
of SDEs, which relies on good strong convergence properties. 

Key words: Stochastic differential equations, non-Lipschitz coefficients, 
Lamperti Transformation, backward Euler-Maruyama scheme. 
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1 Introduction 

The goal of this paper is to derive an efficient numerical approximation for one- 
dimensional SDEs which take values in a domain and have non-Lipschitz drift 
or diffusion coefficients. Typical examples of such SDEs are the Cox-Ingersoll- 
Ross process (CIR), the CEV model, the Wright-Fisher diffusion, where the 
main difficulty is the sublinearity of the diffusion coefficient. Furthermore, the 
approach developed in this paper can be also applied to SDEs with superlinear 
coefficients. Prominent examples are here the Heston 3/2-volatility process 
and the Ait-Shalia model. All the mentioned processes play an important role 
in mathematical finance and bio-mathematical applications. Our key idea is 
to transform the original SDE using the Lamperti transformation into a SDE 
with constant diffusion coefficient, see e.g. [H]. The transformed SDE is then 
approximated by a backward (also called drift-implicit) Euler-Maruyama scheme 
(BEM) and transforming back yields an approximation scheme for the original 
SDE. This strategy was found successful in a recent work [5] for the CIR process, 
where the authors proved that the piecewise linearly interpolated BEM scheme 
strongly converges with rate one half (up to a log-term) with respect to a uniform 
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X p -error criterion. This strategy was also suggested by Alfonsi in pQ. Here, we 
extend that work in several ways: 

• Considering the maximum error in the discretization points, we prove that 
the drift-implicit Euler-Maruyama scheme for the CIR process strongly 
converges with rate one under slightly more restrictive conditions on the 
parameters of the process than in [6] . 

• We provide a general framework for the strong order one convergence of 
the BEM scheme for SDEs with constant diffusion and one-sided Lipschitz 
drift coefficients. 

• Using this framework we present a detailed convergence analysis for several 
SDEs with sub- and super-linear coefficients. 

• We also show that BEM for the transformed SDE is closely related to a 
drift-implicit Milstein scheme for the original SDE, which has been intro- 
duced in [13] . In the case of the CIR process we provide a sharp error 
analysis for this scheme. 

Independently of and simultaneously to the research presented in this paper, 
the same approach was also used by Alfonsi in [5] to derive strong order one 
convergence of the BEM scheme for the CIR and the CEV process. See Remark 
12.91 for a discussion. 

To illustrate the main difficulties and also our main idea let us consider the 
CIR process 



with 2k8 > a 2 . It is a simple implication of the Feller test that the solution of 
equation (|TJ) is strictly positive when 2k9 > a 2 and y(0) > 0. This SDE is often 
used in mathematical finance for interest rate or stochastic volatility models. 
However scalar SDEs with square root diffusion coefficients appear not only in 
the financial literature but belong to the most fundamental SDEs as they are 
an approximation to Markov jump processes [7]. 

Once we attempt to simulate ([1]) using classical discretization methods, see 
e.g. [53], we face two difficulties: 

• In general, these methods do not preserve positivity and therefore are not 
well defined when directly applied to equation JT]); 

• The diffusion term is not Lipschitz continuous and therefore standard 
assumptions required for weak and strong convergence, see e.g. [24] . do 



Consequently, a considerable amount of research was devoted to the numer- 
ical approximation of this equation, see [TTJ [TJ [3TJ |H 13] , to mention a few. 
However no strong convergence of order one results have been obtained so far up 
to best of our knowledge. For a comparison of the different proposed schemes 
based on simulation studies, see [Tl [26], 

Our approach is based on a suitable transformation of the CIR process. 
Applying Ito's formula to x(t) = J y(t) gives 





(1) 



not hold. 





(2) 
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Zhu [3T] pointed out that the drawback of the transformed equation is that the 
new mean level (8 — j^)x(t)~ 1 is stochastic and that a naive Euler discretiza- 
tion cannot capture the erratic behavior of a;(t) _1 -term, although almost sure 
convergence of this method holds true, see [TH]. The weakness of a naive Euler 
discretization is that its transition density is Gaussian and therefore its moments 
explode due to the x(i) _1 -term. 

On the other hand, Alfonsi showed in [TJ that the BEM applied to ([T]) pre- 
serves positivity of the solution and also monotonicity with respect to the initial 
value. Moreover, his simulation studies indicated good convergence properties 
of this scheme. In this paper we follow the recent result by Dereich et al. [6] , 
where it was shown that the piecewise linear interpolation of BEM applied to 
(TT|) , strongly converges with a rate one half (up to a log-term) . 

Given any step size At > 0, the BEM scheme has the form 

X k+1 =X k + -K(6 v X^ 1 -X k+1 )At+~crAw k+1 , fc = 0,l,..., (3) 
X = x(0) 

with 

Aw k+1 = w((k + l)Ai) - w(kAt), k = 0, 1, . . . 

and 




We will establish a strong convergence of order one for the maximum L p -distance 
in the discretization points between (TTJ) and ([3]), see Section 3. For example for 
the L 2 -distance we obtain 

2k9 

E max \Xu - x(kAt)\ 2 < C 2 ■ At 2 for > 3. 

k=0,...,\T/At] a 2 

As a consequence we also obtain the same convergence order for the approxi- 
mation of the original CIR process by X 2 , k = 0, 1, . . . , \T/At] . 

In this paper we will show that the above idea naturally extends to many 
types of SDEs with non-Lipschitz coefficients. The combination of the Lamperti 
transformation and the backward Euler scheme enables us to analyse the L p - 
convergence rates for many scalar SDEs encountered in practice. In particular, 
transforming BEM back we obtain an order one scheme for the original SDE, 
which is close to a Milstein-type scheme, see Section 4. Hence our approach 
turns out to be a new method for deriving numerical methods with strong 
order one convergence for SDEs with non-Lipschitz coefficients. Although strong 
convergence of backward schemes for SDEs with non-Lipschitz coefficients has 
already been analysed in the literature and their convergence for models as the 
Ait-Sahalia and the Heston 3/2- volatility was obtained, see [T31 [57J , schemes 
with strong convergence order one have not been established yet in this setting. 

Another motivation for our work are results by Giles [51 who showed 
that for optimal MLMC simulations one should use discretization schemes with 
strong convergence order one. Note that strong convergence of the discretization 
scheme used for the MLMC simulations seems to be not only a sufficient but 
also a necessary condition |16j . 
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The remainder of this paper is structured as follows. In the next section, 
we provide a general convergence result for the BEM method applied to scalar 
SDEs with additive noise. Section 3 contains the results for our examples, i.e. 
the CIR, CEV, Ait-Sahalia, 3/2-Heston volatility and Wright-Fischer SDEs. In 
Section 4, we provide the relation of BEM and a drift-implicit Milstein scheme 
and give an error analysis for the case of the CIR process. The last section 
contains a short discussion. 



2 The BEM scheme for SDEs with additive noise 
2.1 Preliminaries 

Let D = (I, r), where — oo < I < r < oo, and let a, b : D — > D be continuously 
differentiable functions. Moreover, let (fi, J 7 , (J 7 i)t>o>P) be a filtered probability 
space and w(t), t > 0, a standard (J r t ) t >o-Brownian motion. We begin with the 
SDE 

dy(t)=a(y(t))dt + b(y(t))dw(t), t>0, y(0) € D, (4) 
and assume that it has a unique strong solution with 

P(y(t) &D, t > 0) = 1. 
If b(x) > for all x € D, then we can use the Lamperti-type transformation 

F{x) = x l°W) dy 

for some A > 0. Ito's Lemma with x(t) = F(y(i)) gives the transformed SDE 
dx(t) = f(x(t))dt + Xdw(t), t > 0, ar(0) € F(D), 

with 

where F(D) — (F(l),F(r)). Note that the classical Lamperti transformation 
corresponds to A = 1. This transformation allows to shift non-linearities from 
the diffusion coefficient into the drift coefficient. Then, under appropriate as- 
sumptions on / (respectively a and 6), one can apply the backward Euler scheme 

X k+1 =X k + f(X k+1 )At + XAw k+1 (5) 

and derive a strong convergence rate of order one for the maximum L p -error in 
the discretization points, see Theorem 12.71 in the in Section 2.3. 
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2.2 Backward Euler-Maruyama scheme 

In this section we focus on the numerical approximation of 

dx(t) = f(x(t))dt + adw(t), t>0, x(0) = x a (6) 
by the backward Euler-Maruyama scheme 

X k+1 = X k + f(X k+1 )At + aAw k+1 , k = 0,l,..., X = x . (7) 

We will work under the following assumption on the SDE itself: 

Assumption 2.1. Let — oo < a < (3 < oo and assume that SDE has a 
unique strong solution which takes values in the set (a, (3) C K, i.e. 

P(z(t) € (a, /3), t > 0) = 1. (8) 

For the well-definedness of the drift-implicit Euler-Maruyama scheme we 
need the following assumption on the drift-coefficient: 

Assumption 2.2. The function f : (a, 13) — > (a, (3) is continuous. Moreover, 
there exists a constant if el such that 

(x-y)(f(x)-f(y))<K\x-y\ 2 (9) 

for all x, y £ (a, (3). 

The Feller test, see e.g. Theorem V.5.29 in [22 , gives that condition © is 
equivalent to 

lim v(x) — oo, lim v(x) = oo 

x — >a + x— >/3_ 

with 



2 r f^p'^),,,, 

O- 2 J X0 Jxo P'(Ci) 



and the scale function 



p(x) = J exp ( -2 J ^ du J d£, i£ {a, (3). 



I x 

Note that v can be rewritten as 



= A. r [ 2 exp ( -2 / " Mrfu ] dCidC2- 



u 



IX J Xq y JQi " J 

The condition on v now implies (recall that a < xq < /3) 



and 



lim / exp I 2 / Cl ^ du] dd = oo (10) 



lim / rxp ( I ^-du) dC,i =oo. . I ! ) 
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Now consider (fTUf and assume that \imsnp x ^ a+ f(x) ^ oo However, if this 
would be true, the expression in (|10[) would be finite due to the continuity of /. 
Using a similar argument for ([lip we obtain 

limsup/(x) = oo, liminf f(x) = — oo. (12) 

x — >a+ a;— >p_ 

The drift-implicit Euler scheme is well defined if 

X k+1 - f(X k+1 )At = X k + aAw k+1 

has a unique solution for k = 0, 1, . . .. This is guaranteed by the following result: 

Lemma 2.3. Let Assumvtion \2.1\ and \2.2\ hold and let KAt < 1. Moreover set 

G(x) = x - f(x)At, x£(a,P). 

Then for any eel there exists a unique x G (a,/3) such that G(x) = c. 

Proof. The result follows if we can show that the function G is continuous, 
coercive and strictly monotone on (a, /?) (see |30j ) . However, due to Assumption 
12.21 the function G is continuous on (a,f3). Moreover, since 

(x - y)(G(x) - G(y)) = (x - yf - (x - y)(f(x) - f(y))At 
> (1 - K + At){x - yf > 

by © (with K + = max{0, K}) the required strict monotonicity is obtained. 
Finally, (|12[) and the monotonicity imply that 

liminf G(x) = lim G(x) = — oo 

x— >a+ x— ¥a.+ 

and 

limsupG(a;) = lim G(x) = oo, 

x->/3_ x->/3_ 

so the function G is coercive on (a,/3). □ 
Note that for K < there is no restriction on At. 

For completeness, we state here a well known discrete version of Gronwall's 
Lemma: 

Lemma 2.4. Let At > and let g n ,X n € M, n 6 N, and 77 > 6e given. 
Moreover, assume that 1 — ?yAi > and 1 + A„ > 0, n G N. Then, if a„ G M, 
n G N, satisfies ao = and 

a«+i < a„(l + A„) + i]a n+1 At + g n+1 , n = 0, 1, . . ., 

iften i/w's sequence also satisfies 

^ n—l n—1 

a "" TT^a^E^-^Vi n( 1 + A ')' n = 0,l,.... 

Under the above assumptions we have the following moment bounds for the 
SDE and the BEM scheme: 
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Lemma 2.5. Let T > and let Assumvtion \2.1\ and \2.2\ hold. Then we have 

E sup \x{t)\ q < oo 
te[o,r] 

for all q > 1. // additionally IK At < 77 /or some 77 < 1, then for all q > 1 £/iere 
exzsi constants C q > 0, which are independent of At, such that 

E sup \X h \ q < C q . 

k=0,...,\T/A(\ 

Proof. The first assertion can be shown by a straightforward modification of the 
proof of Lemma 3.2 in [12] . 

For the proof of the second assertion, we will denote constants which are 
independent of At and whose particular value is not important by c regardless 
of their value. Due to (|12p and the continuity of / there exists an x* £ (a, (5) 
with f(x*) = 0, so we can rewrite BEM as 

X k+1 - x* - X k ~x* + (f(X k+1 ) - f(x*))At + aAw k+1 . 

Multiplying with X k+ \ — x* and using the one-sided Lipschitz condition on / 
yield that 

(X k+1 - x*f < {{X k - x*) + aAw k+1 ) (X k+1 - x*) + K+At(X k+1 - x*f . 

where K + — max{0, K}. Moreover, using ah < \{a 2 + b 2 ) and rearranging the 
terms gives 

{X k+1 - x*) 2 < 2K+At(X k+1 - x*) 2 + {(X k - x*) + aAw k+1 f 
and hence, with u k = X k — x* , we have 

u 2 k+1 <u\ + 2K+u 2 k+1 At + 2u k Aw k+ i + a 2 \ Aw t+1 1 2 . 
Note that 1 - 2K + At G (0, 1] and 

sup sup (1 - 2K+At)- k < 00. (13) 

Ate(0,r;/(2if+)) fc=l,..., [T/At] 

So, the above discrete version of Gronwall's Lemma now yields 

k 1 

ul - c+ (i- 2 l+At) fc ^ (1 " 2K+At y^ Aw ^ ( i4 ) 



(1 - 2K+At) k ._ n 



5^(1 - 2K+Aty<T 2 (\Aw 3+l \ 2 - A), 



from which we obtain easily by induction that 

sup Ku 2 < 00. 

k=0,...,\T/A] 



7 



Using this it can be easily checked that the processes , i = 1,2, are square- 
integrable martingales with respect to the filtration TkAt, k — 0, 1, . . .. Hence 
Doob's inequality and straightforward calculations give for any q > 1 that 

E sup \M { k l) \ 2q <c- sup E|u fc | 29 

fc=0,...,[T/Al k=0,...,\T/At] 

and 

E sup \M {2) \ 2q < oo. 

fe=0,...,[T/At] 

So using (fTlf and (|T^|) we obtain 

E sup < c + c ■ sup Ew 2<? 

k=0,...,\T/Ai] k=0,...,\T/At~\ 

for Ai < rj/(2K + ) and the assertion follows now by an induction argument in 
q E N. □ 



2.3 The Main Result 

Here we prove our general theorem on the strong convergence of the numerical 
scheme to the solution of SDE @. 

Assumption 2.6. Let T > and p > 2. We assume that the drift coefficient 
f : (a, j3) —> (a, j3) of SDE ([6]) is twice continuously differ entiable and satisfies 



sup E\f'{x(t))\ p + sup E 
te[o,T] te[o,T] 



< oo. 



Theorem 2.7. Let T > 0, p > 2, n £ (0, 1) and Assumptions \jJj [Ql and\2i 
hold. Then, for 2KAt < n, there exists a constant C' p > ( independent of At ) 
such that 



E 



sup \x(kAt) ~X k \ p 

k=0,...,\T/At] 



< C p ■ At P . 



(15) 



Proof. Recall that we will denote constants which are independent of At and 
whose particular value is not important by c regardless of their value. By 
applying Ito's formula on f(x(t)) we have 

Hfc+l)At 

x((k + l)Ai) = x(kAt) + / f(x((k + l)At))dt (16) 

JkAt 

Ak+l)At 

+ a dw(t)+R k+1 

JkAt 



where 



+ 1 



(fe+l)At r (k+l)At 



a 



(ff)(x(s)) + Y f(x( S )))dsdt (17) 



kAt Jt 

Ak+l)At n(k+l)At 

o-l / f'(x(s))dw(s)dt 

JkAt Jt 
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We can decompose R s as R s = Rs + R^ with 



r sAt r sAt 

RU = - / 

J(s-l)At Jt 

rsAt rsAt 

i?( 2 ) = — (j I / f'(x(u))dw(u)dt 

J{s-l)At Jt 



(ff)(x(u)) + —f"{x{u)) ) dudt 



Using equations and ([7]) we have 
x((k + l)At) - X k+1 = x(kAt) - X k + [f(x((k + l)Ai)) - f(X k+1 )]At + R k+1 
and thus 

(x((k + l)Ai) - X fe+1 - [/(z((fc + l)At)) - /(X fe+1 )]At) 2 = (x(kAt) -X k + R k+1 f. 

We arrive at 

\x((k + l)Ai) - X k+1 \ 2 - \x(kAt) - X k \ 2 

= 2{x((k + l)At) - X k+1 )[f(x((k + l)At)) - f(X k+1 )}At 
- [f(x((k + l)At)) - /(X fc+1 )] 2 Ai 2 + 2(a;(fcAt) - X fe )i? fc+1 + 

Using the one-sided Lipschitz assumption on / we obtain 

(1 - 2K+At) \x((k + l)At) - X k+l \ 2 < \x(kAt) ~X k \ 2 + 2(x{kAt) - X k )R k+1 + R 2 k+1 . 

Let us define e k — x(kAt) — X k and jAt — (1 — 2K + At). Note that jAt & (0, 1] 
and 



sup sup 7 A f < oo. 

Ate(o,7 ) /(2_ff+)) fe=i,...,rr/Ai] 



(18) 



Now, Lemma \2 .41 yields 



A— 1 



k-1 



4 <2j2lAt k ZsRs + l+J2"fA t k R 2 + l- 



(19) 



Since Efi^l-TW] = °> 

we have that 



fc-i 



E^esi?^, fe = 0,...,rT/Atl, 



is a martingale and the Burkholder-Davis-Gundy inequality implies that 



E 



sup 
k=o,...,e 



fc-i 



9/2 



<cE(E7i s *|e s | 2 



I? 



(2) 
s+1 



for any q > 1 and £ = 0, . . . , [T/At] . So using the boundedness of 7a* and ([T 
we have 



E 



sup 

k=0,...,l 



k-1 



Y,^sR 



k„ d(2) 
s+1 



s=0 



9/2 



<cE El 



\s=0 



I? 



(2) 
s+1 



for any q > 1. Using this and Jensen's inequality in (1191) we now arrive at 



E 



sup \e k 

k=0,...,l 



2q 



£-1 



R {2) 



<c\T/At'] q / 2 - 1 J2^\ e s\ q 

£-1 

+ c \T/Aty- i Y / n^\ 

s=0 
l-l 

+ c\T/Aty~ 1 J2 E \R s +i\ 2q 



(20) 



H s+1 



Now, Assumption 12.61 Jensen's inequality and the Burkholder-Davis-Gundy in- 
equality give that 



E 



R {2) 



J~kAt 



<c(i + \f(x kA )r)At 3m ^ 



and 

E 



R {1) 

n k+l 



J~kAt 



< C 1 



(ff)(x(kA)) + Y f"(x(kA)) 



(21) 



At 2m (22) 



for all m < p. Thus, the Cauchy-Schwarz inequality and Assumption 12.61 yield 
that 



d(2) 
H s+1 



E 



|e s | 9 E 



R (2) 
K s+1 



At 



Hence Young's inequality implies 



< 



e-i 



c(l + E|e s | 2 ") 1/2 A^/ 2 



cfT/ At] 9 / 2 - 1 £}E|e s |« q < c^^\e s \ 2q At + cAt 2q 



s=0 



s=0 



Similar we also obtain 



i-l 



\T/At\ q - 1 ^W [ e s \ q R% < cJ2^\e s \ 2q At + cAt 2q . 



Since finally 



c\T/At] q - 1 J2 ]S <\ R s+i\ 2q < cAt2q , 



s=0 



by inserting these three estimates in (|2"0"]l we end up with 



E 



sup \e k 
k=o,...,e 



2g 



< cJ2 E \es\ 2q At + cAt 2q 



s=0 



and Gronwall's Lemma completes the proof. 



□ 



The above result and Lemma 12.51 now give convergence (without a rate) in 
all L 9 -norms: 
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Corollary 2.8. Under the assumptions of Theorem 



lim E 



sup \x(kAt) - X k \ q 

k=0,...,\T/At] 



we have 
= 



(23) 



for all q > 1 . 



Remark 2.9. In J^j, independently of the research in this paper, a similar 
result to Theorem \2. 7| is established for the case D — (a,oo) and drift functions 
f : D — > R, which are twice continuously differentiable and satisfy a monotone 
condition (which is equivalent to our one-sided Lipschitz condition). Using a 
continuous extension of BEM Alfonsi obtains the error bound (|15[) under the 
assumption 

\ P/2 




\f'(x(t))\ dt 




U'f)(x(t)) + -f"{x{t)) 



dt 



< oo 



Propo- 



for p > 1. This result is then applied to the CIR and CEV process, i.e. 
sitions [3~T\ and \3.S\ are obtained. 

Note that due to our bound on the inverse moments on the LBE, see the 
subsection below, we are also able to cover SDEs like the the Heston 3/2-volatility 
and the Ait-Sahlia model. Moreover, since we work under the assumption D = 
(a, /3) we can also treat the Wright- Fisher SDE and similar equations. 



2.4 Boundedness of inverse moments of BEM 

If the drift coefficient has an even more specific structure, see the assumption 
below, we can also control the inverse moments of BEM. For the Heston-3/2 
volatility and also the Ait-Sahalia model this will be helpful later on. 

Assumption 2.10. Let a > and assume that the drift coefficient f : (a,/3) —¥ 
(a, ft) has the structure 

/(*) = Jb+Mz), *e(a,/3) 

where 

\h(x)\ < c 2 ■ (I + \xD xe(a,p) 
for some c\,C2 > and mi,ni2 > 0. 

Under the above assumption BEM can be written as 

X k+1 =X k + {ciX-™ 1 + h(X k+1 ))At + aAw k+1 

and we have 

T^T = An ( Xfe +! - Xk - M^+i)Ai - aAw k+1 ) . (24) 
X k+1 ciAt 

Proceeding as in the proof of Theorem 12.71 we also have 

( , 1 >, = ~tt (x(t k+1 ) - x(t k ) - h(x(t k+1 ))At - aAw k+1 - R k+ i) , (25) 
with R k+ i given by (|T7|). This can be used to derive the following result: 
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Lemma 2.11. Let T > and p > 2. Moreover, let the assumptions of Theorem 
2.7\ and let also A ssumvtion \2.lO\ hold. Then there exists constants C^\c^ > 
such that we have 

sup E\X k \- miP < Cf } I 1+ sup E\x(t)\-' mip ) . 

fc=0,...,[T/At| V t£[0,T] / 



E sup \X k \- m ^<C^ 1+E sup |a;(i)|- miP 



Proof. We only prove the second assertion, the proof of the first assertion is 
similar. Using (|2 % (|25 )| . (TT5|) . ([2"Tjl and (JUJ) we have 

E sup |x(t fc )- mi -X^ mi |f <c(f+E sup |h( I (t i ))-h(Jf t )r 
fc=o,...,rr/At] y fe=o,...,[T/At] y 

Since 

- %)| < \Hx)\ + \h(y)\ < ca(l + |zp 2 + fop"), 
Lemma 12.51 and the triangle inequality finish the proof. □ 

In the next Lemma we establish a general a-priori estimate for uniform 
inverse moments of SDE ([5]) with the drift structure imposed by Assumption 

Lemma 2.12. Let p > 2. Let the assumptions of Lemma \2.5\ hold and in 
addition let the drift of SDE © satisfy Assumvtion \2.lU[ Then there exists a 
constant C p > such that we have 

E sup \x(t)\~ p < C p f 1+ sup E\x{t)\- ip+2 A . 

*e[o,T] y te[o,r] y 

Proof. Let {a„} ngN and {/3n} nC N be such that a„ \ a and /3 n f~ j3, when 
n — > oo. Let us choose no such that a„ < x(0) < f3 no . Then for n > no, we 
define the stopping time r„ = {t > : x(i) ^ (a„, /?„)}. By Ito's lemma we have 



where g(x) = -cip|x|-( p+1+mi ) + <72p{p+1) \x\~^ p+2 \ Observe that for mi > 1, 
lim-r^o g{ x ) — ~°° an d lim^^oo g(x) = 0, hence in that case there exists a c > 
such that 

supg(x) < c. 

x>0 

If mi £ (0, 1] then there exists a c > such that 

i5(^)i<c(i + Nr (p+2) ). 



\x(tAr n )\- p = \x(0)\- p + J (g{x{s))-p\x{s)\-^h{x{s))) l [0 , Tn] (s)ds 
-pa /V(s)r (p+1) l[o,.„](s)d W ( S ), 
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By Assumption 12.101 and Burkholder-Davis-Gundy's inequality we have 



E 



sup |a;(tAr n )| p 
te[o,T] 







<\x(0)\- p + 









■ cE 



1/2 



\x( s )\-^ +1 h [0 , Tn] ( s )d s j 

T 

<\x(0)\- p + I c(l + E\x(s)\-( p+2) ^ ds 



+ cE sup \x(tAT n )\- p I |2;(s)r (p+2) ds 
\te[o,r] 



1/2 



Applying Young's inequality to the last summand of the above inequality now 
yields 



E 



sup |ar(t A r n ) | p 
te[o,T] 



<c 1+ sup E\x{t)\-( p+2 ) 
\ te[o,T] 



and Fatou's Lemma completes the proof. 



□ 



3 Examples 

In this section, we will apply our main result to several examples. To simplify the 
presentation, we will denote the numerical method Yfe = F~ 1 (Xk), k = 0, 1, . . ., 
where Xk, k = 0,1,... is given by ([S]), as Lamperti-backward Euler (LBE) 
approximation of SDE (j4]). Moreover, we will say that this method is p-strongly 
convergent with order one, if 

E sup \y{t k ) - Y k \ p < C p ■ At p . 

k=0,...,\T/At] 

Finally, constants whose particular value is not important will be again denoted 
by c. 



3.1 CIR process 

Recall that the Cox-Ingersoll-Ross process is given by the SDE 

dy{t) = K(9-y(t))dt + (7y/y(t)dw(t), t > 0, 2/(0) > 0. (26) 

If 2k9 > a 2 , then we have D = (0, oo) and Assumption (12.11) holds for (a, (3) — 
(0, oo). Moreover, recall that the transformed SDE using F(y) = ^fy reads as 

dx{t) = f{x(t))dt + -adw(t), t > 0, x{0) = ^/y{$) (27) 

with 

f( x ) = l K f(9-^-jx- l - X Y x>0 
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and the BEM scheme is given by 



X k+l = X k + f(X k+1 )At + -aAw k+1 , k = 0, 1, . . 



(28) 



with Xq = ;r(0). Straightforward calculations give 

( X -y)(f( X )-f(y)) < -l K [x-yf, 



x,y>0, 



so Assumption 12.21 holds with K — —n/2. Observe also that 

1 



f'(x) = -^ K (9 v x- 2 + l), 



and 



(/'/)(*) + y /"(*) = - y (^" 3 -x) + \kB v x-*o*. 



So, for Assumption 12.61 to hold we need 



sup E[x(t)- 3p ] = sup E[y(t)-i p ]. 

0<t<T 0<t<T 



Since 



sup E[y(t) q ] < oo for q > - 



2k0 



0<t<T 



(29) 



see e.g. [BJ, Assumption 12.61 and as a consequence Theorem 12 . 71 hold if p < ~ 
In order to approximate the original CIR process observe that 

(x(t k f - Xl) = (x(t k ) + X tk )(x(t k ) - X tk ). 
Let e > such that p(l + e) < f ■ Then Holder's inequality gives 



E 



sup \x{t k y - xi 

0<fc<|"T/At] 



< E 



sup \x(t k ) + X tk \ p ~ 

0<k<\T/At\ 



E 



sup \x(t k )-X tk \ 

0<k<[T/At] 



p(l+e) 



Using Lemma l2~5l we obtain: 

Proposition 3.1. Let T > 

of the CIR process is p-strongly convergent with order one. 



Proposition 3.1. Let T > and 2 < p < §ff ■ Then, the LBE approximation 



3.2 Numerical Experiment 

Note that the unique solution to (|2"8]l is given by 



■V/. i = 2 ZjT^Af ( ^ fe + 2 lcrAu ' fe+1 + V (^ fe + 2 "°' Au ' fc+1 ) + 
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with 6 V = 6 — j-. Hence implicit schemes not necessarily increase the computa- 
tional complexity with comparison to classical explicit procedures. 
In our numerical experiment, we focus on the L 2 -error at the endpoint T = 1, 
so we let 

e A t = E\x(T)~X T \ 2 ■ 

For our numerical experiment we set 9 — 0.125, k — 2, and a — 0.5. This gives 
= 2 and corresponds to the critical parameters for which Dereich et al. [CT 
established strong convergence of order one half for linearly interpolated BEM 
(|28| with respect to the uniform L 2 -error criteria. Although, theoretical results 
obtained in this paper impose slightly more restrictive conditions for parameters 
than those in Dereich ct al., performed numerical experiment suggests that for 
practical simulations condition 1 < p < 2^ suffices for uniform-L p convergence 
with order one. Although an explicit solution to (|27p is unknown, Theorem 
12.71 guarantees that BEM strongly converges to the true solution. Therefore, 
it is reasonable to take BEM with a very small time step, we choose At — 
2 -15 , as a reference solution. We then compare it to BEM evaluated with 
(2 4 Ai, 2 5 At, 2 6 At, 2 7 At) in order to estimate the rate of the L 2 -convergence, 
where we estimate eAt by a Monte-Carlo procedure, i.e. 

1 1Q3 i 2 

i=l 

Here (T),X^ are iid copies of x{T),Xt- We plot eAt against At on a log-log 
scale, i.e. if we assume that a power law relation eAt = CAt q holds for some 
constant C and q, then we have log eAt = log C + q log At. For our simulation, 
a least squares fit for logC and q yields the value 1.9332 for q with a least 
square residual of 0.016. Hence, our results are consistent with strong order of 
convergence equal to one. 



SVslrong.m 



— * — Backward EM Scheme 
- - Ref. slope At 2 




Figure 1: Strong error plot for backward Euler scheme applied to CIR process. 
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3.3 Heston 3/2-volatility 

In [TU] the inverse of a CIR process is used as a stochastic volatility process, 
which gives the so-called Heston-3/2-volatility 



dy(t) = c x y{t){c 2 - y(t)) dt + c 3 y(i) 3/2 dw(t), t > 0, 
where Ci,C2, C3 > 0. Using F(y) = y^ 1 ^ 2 leads to 



1/(0) > (30) 



dx(t) 



3d 



x{t)- 1 



C\C 2 



x(t) dt 



C3 



dw(t), 



(31) 



which coincides with (|27[) if we use a reflected Brownian motion, i.e. — w, which 
is still a Brownian motion, and 



o = c 3 , 
Hence we have the relation 



1 

(>2 



C\C 2 



K = CiC 2 . 



k6 
T2 



1 



£1 

r 2 ' 



CT c 3 

so Theorem [2~71 holds here forp < |(l + %). Note that the Heston-3/2- volatility 

c 3 

is one of the SDEs which does not have finite moments of any order. As the 
inverse of the CIR process it has finite g-moments up to order q < 2 + ^ . 

c 3 

Now, for transforming back we have to control the inverse moments of the 
BEM scheme for CIR. Here Lemma 12.111 and 12.121 give 



sup 

fc=0,...,[T/At| 



\x k \- p <Cjl 



- sup E\x(t)\-^ +2 '> 
te[o,T] 



for p < |(1 + %). From the analysis of the CIR process we have that 
sup E\x(t)\- (p+2 ) < 00 for p < 4 ( i + % 

t€[0,T] \2 C 3 

To establish the convergence result for the LBM for the Heston-3/2 volatility 
note that 



1 



1 



X k + x{t k ) 

x 2 x{t k y 



x{t k y 

Finally using Holder's inequality gives 



\X k - x(t k )\ < c(\X k \~ 3 + \x{t k )\- 3 )\X k - x{t k )\. 



E 



sup 

fe=0,...,|"T/At] 



< c 



1 

^1 



1 



x{t k ) 2 



3/4 



3/4N 



E sup \X k \~ A P 

k=0,...,\T/At] , 



: sup \x(t k )\~ 

fe=0,...,[T/At] 



x f E sup |X fe -.x(4.)| 4p l 
V fe=o,...,fr/Ai] / 



Hence we obtain: 
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Proposition 3.2. If 1 < p < | + |%, i/ien i/ie approximation of the 

Heston-3/2 process is p-strongly convergent with order one. 

3.4 CEV process 

Another popular SDE in finance is the mean reverting constant elasticity of 
variance process ([5]) given by 

dy{t) = k(9 - y{t))dt + ay{t) a dw(t) (32) 

where 0.5 < a < 1, k, 6, a > 0. By the Feller test we have D = (0, oo) . Applying 
Ito's formula to the function F(y(t)) — y(t) 1_Q we obtain that Assumptions 12. II 
holds with (a,f3) = (0,oo) and 

dx(t) = f{x{t))dt + (1 - a)crdw(t) (33) 

with 

( a QffJ 2 — 1\ 

f[x) = (1 — a) I km !- q - ie , x > 0. 



Again we need to check the remaining assumptions. Since a > 0.5, and conse- 
quently > 2, we have for 

2 

f'(x) = -aKfe _Tr ° — (1 — a)n + (1 — a)— — x~ 2 , x > 

that 

lim f'(x) = — oo, lim f'(x) = —(1 — a)«; 
and hence there exists a c > such that 

sup / (x) < c. 

x>0 

Now the mean value theorem implies 

( x - y)(f( x ) - filj)) <c\x-y\ 2 , x,y > 0, 

i.e. the drift coefficient is one-sided Lipschitz. We have moreover that 

f"(x) = -?—k6x-t^ - (1 - a)aa 2 x~ 3 . 
1 — a 

However, from [4] it is known that 

sup E \y(t)\ p < oo for any p e R, T > 0, (34) 

0<t<T 

and therefore also Assumption 12 .61 holds and Theorem l2.7l can be applied for any 
P > 1. For the back transformation note that the mean value theorem yields 

\( X Th -ytt)\ < (a* +y^)\x- y\, 

1 — a 

Using Lemma [2~5l we have: 

Proposition 3.3. Let p > 1. The LBE approximation of the CEV process is 
p-strongly convergent with order one. 
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3.5 Wright-Fisher Diffusion 

The Wright-Fisher SDE that originated from mathematical biology, [7] and re- 
cently is also gaining popularity in mathematical finance [21 [55] reads as 



dy{t) = (a - by{t))dt + jy/\y{t){l - y(t))\dw(t), t>0, y G (0, 1) 
with a, b, 7 > 0. If 

2a 2(6 -a) 

— > 1 and V _ ; > 1, (35 

then this SDE has a unique strong solution with 

¥(y(t) e (0,1), t>0) = l, 

see [55]. Using 

F(y) = 2arcsin(Vy), 2/ € (0, 1) 

we obtain 

dx(t) = f(x(t))dt + 7 dw(t), t > 0, x = 2 arcsin(Y / yo) 

with 

/(*) = (a - cot (f)-(6-a- tan (|) , x€ (0, tt). 

Since 

A.) + , (f))-5(» — + - '(!)) 

the mean value theorem implies that Assumption 12.21 is satisfied with if = 0. 
Now note that 



/"(*) = i(« -£) cot (f)(l + cot* (f 



1 - a - -t-) tan (^\ (l + tan 2 1 ' 



2 V 4 J V2V V V2 

and 

(/'/)(*(*)) + = (/7)(2arcsin( 2/ (t))) + ^/"(2arcsin( 1 /(t))). 

Since 



cot(arcsin(j/(t))) = y 1 ^-p , tan(arcsin(y(t))) = y - 

and € (0,1) we obtain 

\f\x[t))\<c{l + \y{t)\- l + \l-y{t)\- 1 ) 

and 



(/7)Wi)) + y/"(^W) 



< 



C (l + b(t)|- 3/2 + |l-y(t)|- a/2 ) 
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for some constant c > 0, depending only on a, 6,7 > 0. Using Theorem 4.1 in 
[15] and establishing uniform convergence of the given series expressions in t 
using asymptotic bounds on the Jacobi polynomials we have that 

sup E\y(t)\ qi < 00 
te[o,T] 

if qi > - 2§ and 

sup E\l-y{t)\ q2 < 00 
te[o,T] 

2(fc-q) 

4 



if q2 > — ^2 ■ So Assumption 12.61 is satisfied if 



3", 



min{a, b — a} > p. 



Now Theorem 12.71 gives 



E 



sup \x(t k )-X k \ 

k=0,...,fT/At] 



< C p ■ At p 



for p < — - minia, b — a} 



and transforming back yields 



E 



sup 

fe=0,...,fT/Atl 



y(t k ) - sin 2 (-^j 



< C p ■ At p . 



under the same assumption on the parameters since the sin-function is bounded. 



Proposition 3.4. Let 2 < p < min{a, b — a}. Then the LBE approximatic 
of the Wright-Fisher process is p-strongly convergent with order one. 



3.6 Ait-Sahalia model 

Higham et al. analysed in [55] a backward Euler scheme for the Ait-Sahalia 
interest rate model 

dy(t) = (a-njity 1 -a + a iy (t) - a 2 y(t) r )dt + ay(t) p dw(t), (36) 

where a_i, ao, a, a±, a.2, o~ are positive constants and p, r > 1. In [29) it was 
established that 

p(»(t) e (0,00), t>o) = i. 

Under the assumption r+1 > 2p Higham et al. proved uniform L p -convergence 
for any p > 2 the backward Euler method directly applied to (j36|) . However, 
their results did not reveal a rate of convergence. Here, using the Lamperti 
transformation approach we construct a scheme that strongly converges with 
rate one. We focus on the critical case with r = 2 and p = 1.5 which was not 
covered in [29 . Using F(y) = y^ 1 ^ 2 we obtain 

dxif) — f(x)dt — -adw(t). 
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with 



f(x) = Q«2 + jjcr 2 ^ x(t) 1 - ^aix(t) + ^a x{t) 3 - ^a- 1 x(t) 5 , x > 0. 



We have for 



that 



/'(•'■) = ( \ol 2 + ^cr 2 j x(t) 2 - iai + ^a x(i) 2 - |a-i:r(i) 4 , 



lim /'(x) = — oo, lira f (x) = — oo 

x— >0 2:— >oo ' 



and hence there exists a c > such that 

sup f'(x) < c. 

x>0 ' 

Now the mean value theorem implies 

{x-y){f(x) - f(y)) <c\x-y\ , x,y>0, 
i.e. the drift coefhcient is one-sided Lipschitz. We have moreover that 

f"(x) = (a 2 + \A x(t)- 3 + 3a x(t) - 10a_ix(t) 3 , 



<c(l + \x(t)\-* + \x(t)\ m ). 



and 

(f'f)(x(t)) + ^f"(x(t)) 
Straightforward computations also give that 

f(x)>g(x), 16(0,21(0)] 

where 



9(x) ■■= ( ^ct 2 + ^<r 2 ) x 1 - ^-a 2 (3x 



with 



= h 16 x(0) . 

a 2 a 2 



Now a comparison result for SDEs, see e.g. Proposition V.2.18 and Exercise 
V.2.19 in [22], yields that almost surely 



where 



and 



x(t)>x^(t), (e[0,rW), 



dx£\t) = g{x [1) {t))dt - -adw(t), t > 0, x (1) (0) = x(0) 



r (1) = inf{< G [0,T] : x(t) > 2x(0)}. 
Let us define a sequence of stopping times 

= w£{t G [rP" -1 ) , T] : x(t) < 2x(0)}, i = 1, 2, . . . . 
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r (2i+i) = inf { t g [ T (2i) ; T ] . ^ > 2x (0)}, i = 1,2, ... 
and an associated sequence of SDEs 

da;W( t ) = g( x ®(t))dt - ^adw{t), t > t^ 1 \ x {i) (t^ ) = x{t^). 
Using the comparison result for SDEs again, we have that almost surely 
x{t)>x^{t), te[r( 2 '),r (2,+1) ). 



Therefore, using the CIR process as a lower bound, Assumption 12.61 and conse- 
quently Theorem 12 . 71 holds for p < | (l + ^JJ. Proceeding as for the 3/2-model 
we have: 

Proposition 3.5. Let l<p<4 + |^f. The LBE approximation of the Ait- 
Sahalia process with r = 2 and p = 1.5 is p-strongly convergent with order one. 

In the case r + 1 > 2p we know from [33] that 

sup E|x(t)|~ p < oo for all p > 1. 

te[o,T] 

Moreover, the drift coefficient of the transformed SDE behaves at zero like the 
one of a transformed CEV process. A by now standard analysis gives: 

Proposition 3.6. Letp > 1. The LBE approximation of the Ait-Sahalia process 
with r + 1 > 2p is p-strongly convergent with order one. 



4 A Milstein-type scheme for CIR 

In this section we establish a connection between the Lamperti-backward Eulcr 
and a drift-implicit Milstein scheme for the CIR process. We will show that the 
order of convergence of the LBE carries over to a drift-implicit Milstein scheme, 
which has been proposed in [20] and [13] . While strong convergence was shown 
in [23], sharp convergence rates have not been established so far. 
Recall that BEM for the transformed CIR process reads as 

X k+1 = X k + f(X k+1 )At + ^aAw k+1 , fc = 0,l,... 

with 




x > 0. 



Squaring yields the LBE, i.e. 

X 2 k+ i - Xl + k(9 - X 2 k+1 )At + aX k Aw k+1 

+ ^((Aw k+1 ) 2 - At) - (f(X k+1 )) 2 At 2 . 
On the other hand the drift-implicit Milstein scheme for CIR is given by 

2 

Z k+1 = Z k + k{0 - Z k+1 )At + a^fz~ k Aw k+1 + y((Aw H1 ) 2 - At), (37) 
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hence both schemes coincide up to a term of order At 2 . The numerical flows of 
the LBE and Milstein scheme are given by 

^e{x, k) = 1 fx + n9At + cry/x~Aw k+1 + a — ((Aw fc+ i) 2 - 

1 + kAz \ 4 ' 

and 

Mx > k) = TT^At 

It is clear then that 



2 \ 

nOAt + a^Aw k+1 + ^— ((Aw fe+ i) 2 - At) J 



(f>M{x,k) > (f> E (x,k) 
for all x > 0, k = 0, 1, . . .. From [T] we know on the other hand 

4>e(x, k) > 4>E{y, k) for x > y. 

Hence we conclude 

Z k >Xl fc = 0,l,..., (38) 

so the drift-implicit Milstein scheme dominates the Lamperti-Euler method and 
thus preserves positivity. 

To establish the order of ^-convergence for the drift-implicit Milstein scheme 
it is enough to control the difference between the Lamperti-Euler method and 

Lemma 4.1. Let ^ > 3/2. Then there exists a constant c > such that 

sup ElZfc- Xl\ <cAt. (39) 

fc=0,...,["T/At] 

Proof. Let e k = Z k — X 2 and note that e k > by (|3"51) . We have 

e fc+ i = e k - Ke k+1 At + o-{yfz k ' - X k )Aw k+1 + (f(X k+1 )) 2 At 2 . 
Exploiting the independence of X k , Z k and Aro^i it follows 

Eefe +! = T-^a7 (Ee fc +E(/(X fc+1 )) 2 At 2 ) 
1 + nAt v ' 



and consequently 



[T/Atl-i 

sup Ee fc < At ^ E(/(X fe+1 )) 2 At. 



fc=o,...,rr/Ati fe=0 
Due to our assumptions Lemma [2.111 gives that 

sup EX^ 2 < c, 

k=l,...,lT/At] 

which together with 

^c^l + lxl + H" 1 ) 
and Lemma [2~51 shows the assertion. 



□ 
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Using this result we have: 

Proposition 4.2. (i) Let ^ > 3/2. Then, there exists a constant C > such 
that 



sup E\y(kAt) - Z k \ < C ■ At. 

k=0,...,[T/At] 

(ii) Let ^ > 3/2. Then, there exists a constant C > such that 
E sup \y(kAt) ~ Z k \ 2 < C ■ At. 

fc=0,...,|"T/At] 



(40) 



(41) 



Proof, (i) This follows from the triangle inequality and Proposition 13. II 

(ii) Using (|40]) the second assertion can be shown along the lines of the proof 
of Proposition 5.3 in [23., where strong convergence of the drift-implicit Milstein 
scheme (without a convergence rate) was shown. Proceeding as in the proof of 
Proposition 5.3 in [53] we have 

|T/At"|-1 

E sup \y{kAt) - Z k \ 2 <c ( 1 + KAt) 2l E W^t) - Z e \ At 



k=0,...,\T/At~\ 



c E sup 

fe=l,...,[T/At] 



A— 1 



with 

Pk+l = ~K 

So, (gDJ gives 



(fc+l)At 



kAt 



(y(s)-y((k+l)At))ds+a 



(fe+l)At 



(Vy(^-Vy(kMj)dw(s). 



kAt 



E sup \y(kAt) - Z k \ 2 < cAt + cE sup 

fc=0,...,|"T/Atl k=l,...,\T/At\ 



fs—1 1 

^ (1 + nAtf 



-Pt+i 



For the second term straightforward computations yield that 

2 



E sup 

k=l,...,\T/At\ 



fc-1 . 

^ (1 + «At) fc -< 



Pi+i 



< cAt, 



which completes the proof of the proposition. 



□ 



Using the (suboptimal) second estimate of the above Proposition and Lemma 
3.5 of [B] we also obtain a sharp error estimate for the piecewise linear interpo- 
lation of the drift-implicit Milstein scheme, i.e. 

Zt = — ^ -Z k H jj-^Zfc+i, t 6 [t kl t k+ i], 



A 



A 



in a combined L 2 -\\ ■ lira non 



23 



Proposition 4.3. Let ^ > 3/2. Then, there exists a constant C > such that 

_ \ 1/2 

E max \y(t) - Z t \ 2 ) < C ■ *Jlog(At)\ ■ VM, 

f£[0,T] / 

for all A e (0,1/2]. 

The above relation between a drift-implicit Milstein scheme applied to the 
original SDE 

dy(t) = a(y(t))dt + b(y(t))dw(t) 
and the BEM applied to the transformed SDE 

dx(t) = f(x(t))dt + \dw(t) 

with 
and 

F{x) = x fw) dv 

is in fact a particular case of a more general relation. Expanding LBE yields 
F-^Xk+x) = F-^Xk) + hb(F- 1 (X k ))(X h + 1 - X k ) 

+ ^b'(F- 1 (X k ))b(F- 1 (X k ))(X k+1 - X k ) 2 + ..., 

since 

±F-\ X) = { b{ F-\ x)) , jgpF-H*) = ±«{F-^WF-H*)). 
Setting Y k — F^ 1 (X k ) and using ([5]) we have 
Y k+1 =Y k + (a{Y k+1 ) - h'b(Y k+1 )] At + b(Y k )Aw k+1 + \b'b{Y k )\ Aw k+1 \ 2 



+ R k 

with 

Rk = (b(Y k ) - b(Y k+1 ))f(X k+1 )At + j{b'b)(Y k )Aw k+1 f(X k+1 )At 

+ ^(b'b)(Y k )f 2 (X k+1 )At 2 + ... 

So dropping R k and the other higher order terms, we end up with 

Z k+1 = Z k + (a(Z k+1 ) - ~b'b{Z k+1 )j At + b{Z k )Aw k+1 + ~b'b(Z k )\ Aw k+1 \ 2 . 

In the case of (a, (3) — (0, oo) conditions for the well-definedness, stability and 
strong convergence of this scheme are given in |13j . However, the convergence 
rate analysis for the CIR process, where we can exploit (among other things) 
that F^ 1 (x) — x 2 and also the domination property ([38]) . seems not to carry 
over to the general case. 
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5 Conclusion and Discussion 



In this paper we presented a Lamperti-Euler scheme for scalar SDEs which 
take values in a domain D = (l,r) and have non-Lipschitz drift or diffusion 
coefficients. Our strategy is to first use the Lamperti transformation x(t) = 
F(y(t)) (provided that the diffusion coefficient of the original SDE is strictly 
positive on D) and then to approximate the transformed process x(t), < t < T, 
with the backward Eulcr scheme. Transforming back with the inverse Lamperti 
transformation gives an approximation scheme for the original SDE. We also 
pointed out a relation of this scheme to a drift-implicit Milstein scheme. 

The advantages of this Lamperti-Euler method are 

• that it preserves the domain of the original SDE 

• and an available framework which allows to establish strong convergence 
order one for this scheme. 

In particular, we use this framework to obtain such strong convergence re- 
sults for several SDEs with non-Lipschitz coefficients from both mathematical 
finance and bio-mathematics. 

Whether the implicitness of this scheme (which e.g. for the CEV process 
requires solving a non-linear equation) can be avoided by using a tamed Euler 
scheme (as in jTT] for the case D = M.) remains an open question. Open ques- 
tions are also, whether Assumption 12.21 and 12.61 can be formulated in terms of 
conditions on the original coefficients of the SDE, and whether the convergence 
rate for the Lamperti-Euler also carries over in general (and not only for the 
CIR process) to a drift-implicit Milstein scheme. In particular, the last point 
leads to the general question that given a certain numerical approximation (for 
an SDE with non-Lipschitz coefficients) which perturbations do not change its 
convergence properties and also its qualitative properties? We will pursue all 
these topics in our future research. 
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